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The problem of accommodating unknown sensor bias is considered in a direct model 
reference adaptive control (MRAC) setting for state tracking using state feedback. Sensor 
faults can occur during operation, and if the biased state measurements are directly used 
with a standard MRAC control law, neither closed-loop signal boundedness, nor asymp- 
totic tracking can be guaranteed and the resulting tracking errors may be unbounded or 
unacceptably large. A modified MRAC law is proposed, which combines a bias estimator 
with control gain adaptation, and it is shown that signal boundedness can be accomplished, 
although the tracking error may not go to zero. Further, for the case wherein an asymp- 
totically stable sensor bias estimator is available, an MRAC control law is proposed to 
accomplish asymptotic tracking and signal boundedness. Such a sensor bias estimator can 
be designed if additional sensor measurements are available, as illustrated for the case 
wherein bias is present in the rate gyro and airspeed measurements. Numerical example 
results are presented to illustrate each of the schemes. 

Keywords: Sensor fault accommodation, adaptive control, sensor bias. 

I. Introduction 

Control of uncertain systems in the presence of actuator and sensor faults is an important and challenging 
problem. Two traditional approaches to this problem are fault detection, identification (FDI)-based control 
reconfiguration, and direct adaptive control. Although considerable research effort has focused on actuator 
fault accommodation in the past, literature addressing sensor fault accommodation appears to be rather 
limited. However, sensor faults have been determined to be a cause in several aircraft accidents (e.g., radio 
altimeter fault , 1 angle-of-attack sensor fault , 2 airspeed sensor fault 3 ). 

Direct model reference adaptive control (MRAC) methods aim to control uncertain systems using con- 
troller gains that are adaptively adjusted to achieve a performance close to that of a reference model while 
maintaining system stability and close tracking of the reference model response. Direct MRAC has been 
known to be an effective method for state- or output- tracking. Direct MRAC schemes that use state feed- 
back for state-tracking (SFST) have the advantage of simplicity of implementation while achieving effective 
state tracking in the presence of parameter uncertainties . 4 Such schemes have also been extended to the case 
with actuator failures , 5 , 6 . 7 However, aside from actuator faults, sensor faults may corrupt the measurements. 
Any sensor bias prior to the operation can be compensated for by appropriate off-line calibration methods. 
However, unknown biases can appear during operation in sensors such as rate gyros, accelerometers, altime- 
ter, etc. The most common approach to deal with such sensor faults is to use redundant sensor packages. 
However, common-mode failures can occur across all the sensors, and moreover, each sensor may develop a 
different unknown bias. If used directly in a MRAC control law, such offsets in sensor measurements can 
have detrimental effects on closed-loop stability. 

Sensor bias estimation and accommodation has been a topic of considerable research. For example, bias 
estimation for multisensor systems has been investigated in , 8 , 9 . 10 The problem of setpoint tracking with 
sensor bias was discussed in , 11 and MRAC schemes for output feedback for output tracking in the presence 
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of sensor faults have been developed (e.g.,, 1213 ). However, accommodating sensor bias in SFST-MRAC has 
not been investigated to the knowledge of the authors. 

In this paper, it is assumed that all the a priori bias has been compensated for by off-line calibration 
methods, and trim values of the system states are available at the start of the control operation. Then, 
direct MRAC is considered in the presence of any constant unknown sensor bias that may appear during 
the operation. If biased sensor signals are used directly in the MR AC law, signal boundedness and tracking 
cannot be proved, as discussed in Section II. A logical extension of the standard MRAC is proposed in 
Section III, which incorporates sensor bias estimation in the MRAC law, and shows that signal boundedness 
and bounded tracking error can be accomplished. To compensate for the effect of bias, an additional adaptive 
gain is included in this control law. The case wherein a sensor bias estimator is available is considered in 
Section IV. If the sensor bias estimator has asymptotic convergence properties, a modified MRAC law can 
be developed, which can guarantee signal boundedness and asymptotic tracking. As an example, Section 
V considers a longitudinal aircraft dynamics model with constant biases in the pitch-rate measurement 
gyroscope and in the airspeed measurement. Based on the availability of reliable pitch angle and position 
measurements, non-model-based observers are constructed and used in conjunction with adaptive control. 
The resulting MRAC scheme is shown to provide asymptotic tracking and signal boundedness. Section VI 
contains application examples of the proposed schemes, and concluding remarks are given in Section VII. 

II. Sensor Bias and Standard MRAC 

Consider a linear time-invariant plant described by: 

x ( t ) = Ax ( t ) + Bu ( t ) (1) 

y (t) = x(t) + (3 

where A £ R n x " is the system matrix assumed to be unknown, B £ R " x rn is a known input matrix, 
x(t) £ R n is the system state, and u (t) £ R m is the control input, y (t) £ R n is the available state 
measurement with an unknown constant bias (3 £ R n . 

The objective is to design an adaptive feedback control law using the available measurement y (f) with 
unknown bias (3, such that closed-loop signal boundedness is ensured and the system state x ( t ) tracks the 
state of a reference model described by 

fm (t) = A m x m (f) + B m r (t) (2) 

where x m € R n is the reference model state, A m £ R nxn , B m £ R nxm , and r(t) £ R mr (1 < m r < m) is a 
bounded reference input used in system operation (e.g., pilot input in the case of aircraft). 

It is assumed that the system (A, R) and the reference model ( A mi B m ) satisfy the SFST matching 
conditions, i.e., there exist gains K\ £ R nxm , K 2 £ R mxm r such that 

A m =A+ BKj ; B m = BK 2 . (3) 

The reference model is customarily designed to capture the desired closed-loop response of the plant. For 
example, the reference model may be designed using optimal and robust control methods such as LQR, U 2 , 
or Hoc methods. For the adaptive control scheme, only A m and B need to be known. Moreover, there exist 
positive definite matrices P = P T , Q = Q T £ R nxn , such that the following Lyapunov inequality is satisfied: 

A^P + PA m < -Q. (4) 

The standard MRAC SFST control law is given by 

u = Ki y + K 2 r (5) 

where K\ ( t ) £ R nxm , K 2 ( t ) € R mxm are the adaptive gains. Substituting (5) into (1), the closed-loop 
plant is given by: 

x = Ax + B (Ki y + K 2 r^ 

= Ax + B (Kf y + K 2 r) + B (f<^y + Ayr) 

= (A + RATf ) x + B y + AT 2 r) + BK 2 r + BI< f (3 (6) 

2 of 18 


American Institute of Aeronautics and Astronautics 



where the fact that y = x + 0 was used, and K\ = K\ — K\\ K 2 = I\ 2 — ■ Using the matching conditions 

(3) in (6), we get 

x = A m x + B m r + B (kj y + K 2 r ) + Bl<[ 0. (7) 

The last term in (7) represents a bias disturbance term, which will need to be compensated for by the 
adaptive control law. Defining the tracking error e = x — x m , and using (2) and (7), we have 

e = A m e + B (RJ y + K 2 r ) + BK? (3. 

Due the presence of unknown sensor bias, x(t) cannot be measured; therefore the tracking error e (t) is 
not available. Thus, the standard MRAC control law in (5) cannot guarantee signal boundedness (including 
boundedness of the tracking error) in the presence of sensor bias in state feedback. The next section proposes 
a modified adaptive control law that aims to estimate and compensate for the sensor bias. 

III. Modified MRAC with Sensor Bias Estimation 

In this section, a bias estimator is proposed as a part of the adaptive control law that uses the available 
(biased) measurements, and can provide bounded state tracking error along with signal boundedness. 

Let 0 (t) denote an estimate of the unknown sensor bias (3. Using (3 , define the ‘corrected’ state x(t) £ R n 
as 

x = y-0. (8) 

Therefore, 

x = x + (3 — 1 3 = x + 0. (9) 

where (3 = (3 — (3. Design an adaptive control law as 

u = Ki y + K 2 r + k 3 (10) 

where K\ (t) € R nxm , K 2 (t) £ R rn x mr , and k 3 (t) £ R m are the adaptive gains. Therefore, the closed-loop 
corrected-state equation is 

x = Ax + B [kj y + K 2 r + k 3 S j + P 

= Ax + B (Ki y + K 2 r + k 3 ) + B (k^y + k 2 r + k^j + (3 

= (A + BK f) x + B ( K l y + K 2 r + fc 3 ) + BK 2 r + BK f 0 + Bk 3 + 0 (11) 

where K\ = K\ — K\, K 2 = K 2 — K 2 , and k 3 = k 3 — k 3 . The following matching conditions are assumed to 

be satisfied for some ideal gains K\ , K 2 , k 3 : 

A m = A + BKj, B m = BK 2 , BKl0 = -Bk 3 . (12) 

Using (12) and (9) in (11), we get 

i = A m x + B m r + B (k^ y + k 2 r + k^j - A m 0 + 0. (13) 

Define a measurable auxiliary error signal e(t) £ R n as 

e = x-x m . (14) 

Therefore, from (9), we have 

e = x — x m + 0 = e + 0. (15) 

Differentiating (14) with respect to time, the closed-loop auxiliary error system can be expressed as 

e = x — x m . (16) 
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Substituting (13) and (2) into (16) yields 

e = A m e + B ^ it iU + k 2 r + kkj — A m (3 + /3. (17) 

The following theorem gives adaptive gain update and bias estimation laws that guarantee closed-loop signal 
boundedness as well as bounded tracking error. 

Theorem 1: For the system given by (1), (2); the adaptive controller (10), the gain adaptation laws 

K\ = -r iye T PB 

K 2 = -Y 2 B T Per T (18) 

k = -r 3 B T Pe 

where Tj G R nxn , T 2 G R mxm , T 3 G R rn x m are constant symmetric positive definite matrices and P was 
defined in (4), and the bias estimation law 

b = -vP^AlPe (19) 

where rj € R is a tunable positive constant gain, guarantee that all the closed-loop signals including the 
adaptive gains and bias estimate are bounded and the tracking error e (t) is bounded. 

Proof: Define 

n m r 1 

V = e T Pe + Y, KlV^k u + £ K^k^ + kj P~P (20) 

i i ^ 

where the subscript i denotes the i th column of K 1 , A' 2 . Differentiating (20) with respect to time, using (4), 
(17), properties of matrix trace, and the gain update laws in (18), the following expression is obtained upon 
simplification: 

V < —e T Qe - 2e T PA. m (3 - 2e T Pf3 - -f3 T P/3. (21) 

V 

Using the bias estimation law of (19) in (21), we get 

V < —e T Qe - 2e T PA m p + 2 V e T A^Pe + 2(3 T Aj n Pe 
= — e T Qe — rje T Qe = — (1 + rj) e T Qe. 

Therefore, V < 0, i.e. , V(T) is bounded for all T, and e(t), /3(f) , y(t), k\, k 2 , k are all bounded and 
e(t), (3 (t) G L 2 . From (17), (19) and closed-loop signal boundedness, we have e (t) , (3 G L°°, therefore 
limt^oo e(f) = 0. Because e(t) G L 2 n L°° and y(t) G L°°, K u , K 2 i , ^3 G L 2 n L°° . It can be verified 

that Ku, k 2i , k 3 G L°°, therefore, Hindoo I\u(t) = 0, K 2i (t ) = 0, lim^oo fc 3 ( f) = 0. That is, all 

signals and estimates are bounded, and lim t ^ 00 (a: — x m ) = 0 . ■ 

If sufficient persistent excitation is present, (3 (t) should approach 0 as t — > 00 , and therefore e (f) — > 0 as 
t — > 00 . 

The adaptive control law in Theorem 1 guarantees stability (signal boundedness) and bounded tracking 
error. However, it would be desirable to achieve asymptotic tracking. In an effort to accomplish asymptotic 
tracking, the next section addresses the case when a separate asymptotically stable bias estimator is available 

(e-g-, [14]). 


IV. MRAC with Asymptotic Bias Estimation 


Suppose an asymptotically stable bias estimator is available. This section illustrates how such an esti- 
mator can be employed in a MRAC setting to achieve asymptotic state tracking. 

Consider a bias estimator with estimation error dynamics of the form 


where (3 (f) is an estimate of (3, and 


(3 = Apfi 


P(t)=f3-P(t) 


(22) 
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is the estimation error. Ap £ Rn><n j s a known asymptotically stable matrix, such that lim^oo/? (t) = 0. 
Defining the adaptive control law as in (10) and proceeding as in Section III, we obtain (17). The following 
result gives an adaptive control law that guarantees signal boundedness and asymptotic tracking. 
Theorem 2: For the system given by (1), (2); the adaptive controller (10), the bias estimator (22), and the 
gain adaptation laws 


ki = -r iye T PB 

K 2 = -Y 2 B T Per T (23) 

k = -r 3 B T Pe 

where Ti £ R nxn , T 2 £ R mxrn , T 3 £ R mxm are constant symmetric positive definite matrices and P 
was defined in (4), guarantee that all the closed- loop signals including adaptive gains are bounded and the 
tracking error e (t) — > 0 as t — » oo. 

Proof: Define 

n m r 

V = e T Pe + J2 KlV^k u + ^ K 2i + kj T^k 3 + 0 T Pp[3 (24) 

i i 

where the subscript i denotes the i th column of Ri, I\ 2 \ and Pp = Pj £ R nxn is a positive definite solution 
to the Lyapunov inequality 

A T pPp + PpAp<-Qp (25) 

for some positive definite matrix Qp = Q J £ R nxn . Differentiating (24) with respect to time, using (4), 
(17), (22), (25), and properties of matrix trace, the following expression is obtained upon simplification: 


V < -e T Qe-2e T P{A m - Ap)0- (3 T QpP + 2k^B T Pe + TJ 1 k 3 } 


+2 Tr 


kf { ye T PB + T 




2 Tr 


RJ { B T Per T + 


Using the gain update laws (23) in (26), we get 


V < -e T Qe - 2e T P (A m - Ap) 0 - 0 T Qp0 

= —z T Qz 


where z £ R 2n is defined as 
and Q £ R 2nx2n is defined as 


lT 


e (3 


Q = 


Q 


P (A m — Ap) 


[(Am — Ap) P Qp 
Since Q is positive definite, the matrix Q is positive definite iff 

Qp - (Am - Apf PQ~ l P (Am - Ap) > 0. 


' R‘l } 


(26) 


(27) 


Qp can be chosen such that (27) is satisfied; therefore, V < 0, i.e., V(T) is bounded for all T, and 
e(t), $(t) , y(t), R\ , R 2 , k 3 are all bounded and e(t), (3(t) £ L 2 . From (17), (22) and closed-loop signal 
boundedness, we have e(t) , 0 £ L°°, therefore lim^oo e(t) = 0. From (22), lim^oo/3 (t) = 0. Therefore, 
from (15), lim t ^o 0 e(f) = 0. Because e(t) £ L 2 C I L°° and y(t ) £ L °° , Ku, K 2i , k 3 £ L 2 C\ L°°. It can 

be verified that I<u, R 2i , k 3 £ L°°, therefore, lim^oo Ru (t) = 0, Hindoo I\ 2i (t) = 0, linp^oo k 3 (t) = 0. 
That is, all signals and estimates are bounded, and lim^o^a: — x m ) = 0. ■ 
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V. Non-Model-Based Observers for Longitudinal Aircraft Dynamics 


In this section, we illustrate how sensor bias estimators can be developed and used in conjunction with 
MRAC. The bias estimates are an outcome of non-model-based observers, which additionally, also generate 
estimates of other measurable states. It is shown, how these additional state estimates, and the bias estimates 
are combined with MRAC to yield asymptotic state tracking even in the presence of unknown sensor bias. 
Specifically, we consider a fourth-order longitudinal dynamic model of an aircraft. The state variables 
are: true airspeed v(t), angle-of-attack a(t), pitch angle 6(t), and pitch rate q(t). Thus, the state is x = 
r R 


v a 9 q 


V.A. Rate Gyro Bias Observer 

The pitch rate is usually measured by a rate gyro which is typically subject to bias buildup. Differentiating 
pitch angle to obtain the pitch rate is not advisable due to the presence of noise. Consider that the available 
state measurement (y = x + (3) has an unknown constant bias in the pitch rate measurement, such that 


0 


0 0 0 



(28) 


where b is an unknown constant. Thus, the measured pitch rate ( q rn (t)) can be expressed as 

q m = q + b = d + b. 

Therefore, the actual pitch rate (q = 6) is 

9 = —6 + q m . (29) 

The attitude (pitch angle) sensor usually consists of a two degrees-of-freedom (DOF) vertical gyro. The 
measurements degrade when the gimbal axes are not orthogonal, or a gimbal lock occurs when the spin 
axis is aligned with the outer gimbal axis. Therefore, it may be preferable to use other methods such as 
differential GPS signals to measure the aircraft attitude. Suppose an independent bias-free measurement 
y a (t) of the pitch angle is available as 

Vo = 9. 

Since, the unknown bias b is constant, the Eq. in (29) can be augmented with the equation: 6=0. This 
yields 


±o — A 0 x 0 T B 0 q m 
yo — C 0 x 0 


where 


e 

5 A 0 — 

0 -1 


1 

; C„ = 


; B 0 = 

1 0 

b 


0 0 


0 




Since (C 0 ,A 0 ) is observable, an estimator can be constructed as 

x 0 = A 0 x 0 + B 0 q m + H (y 0 - y Q ) 
where H £ R 2 is an observer gain (e.g., Kalman gain) 


H = 


hi 



Using (30) and (32), we get the following estimation error dynamics 

x 0 (A 0 HC 0 ) x 0 

where x a (t) = x a (t) — x Q ( t ) is the estimation error. Using (31) and (33) in (34), we have 


1*1 


1 

1 

-1 


1 

•1-0 1 


1 

1 

to 

1 

o 


1 

i 


(30) 

(31) 

(32) 

(33) 

(34) 

(35) 
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Since, the actual 6 is measurable, the initial estimate of 9 (0) can be initialized to its true value, such that 
9 (0) = 0, while b (0) can be chosen arbitrarily. 

The control law is defined as in (10), such that using the matching conditions in (12), the following 
closed-loop corrected-state equation is obtained: 

x = A m x + B m r + B [kj y + k 2 r + k 3 ) - A m (3 + (3 (36) 


where 


/? = 


1 T 


0 0 0 b 


Substituting from (2) and (36) into (16), the following closed-loop error system is obtained: 


(37) 


e — A m e + B ^ K i y + k 2 r + k^ — A m f3 + f3. 
Using (37) with (35), (38) can be expressed as 

e = A. m e + B (k^ y + k 2 r + k 3^ — A m k> — h 2 9 


(38) 


(39) 


where A m 4 represents the 4 th column of A m , and h 2 = 
such that 


0 0 0 h 2 


. Define a constant matrix M £ i? 4x2 


M = 


h 2 A. m 4 


Then (39) can be expressed as 


e = A m e + B ^ k"[y+ I< 2 r + - Mx 0 . 


(40) 


The following result shows signal boundedness and asymptotic tracking using such a bias estimator. 
Theorem 3: For the system given by (1), (2); the adaptive controller (10), the gain adaptation laws in (18), 
and the non-model-based observer in (32) guarantee that all the closed-loop signals including the adaptive 
gains are bounded and the tracking error e (t) — * 0 as t — > 00. 

Proof: Define 


It lll r 

V = e T Pe + Y, Klk^ku + Y Kl T 2 X k 2l + k%T 3 % + x T 0 P 0 x 0 (41) 

i i 

where the subscript i denotes the i th column of K\, k 2 \ and P 0 = Pj £ i? 2x2 is a positive definite solution 
to the Lyapunov inequality 

A t 0 P 0 + P 0 A 0 <-Q 0 (42) 

for some positive definite matrix Q a = Q ° £ R 2x2 . Differentiating (41) with respect to time, using (4), (40), 
(42), properties of matrix trace, and the gain update laws in (23), the following expression is obtained upon 
simplification: 


V < —e T Qe — 2e T PMx 0 — XgQ 0 x 0 
= —z t Qz 


where z £ R 6 is defined as z 


e 



and Q £ R 6x6 is defined as 


Q PM 
(• PMf Q 0 


Since Q is positive definite, therefore V < 0 iff Q a — [PM) T Q~ l PM is positive definite. Q$ can be chosen 
such that Q is positive definite. The rest of the proof is similar to Theorem 2. Therefore, all signals and 
estimates are bounded, and lim t ^ 00 (a; — x rn ) = 0. ■ 
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V.B. Velocity Bias Observer 

Similar to the non-model-based observer designed to estimate the bias in pitch rate, an observer can be built 
to estimate the bias in true airspeed if reliable ground position and wind speed measurements are available. 
Consider that the available state measurement (y = x + /3) has an unknown constant bias in the true airspeed 
measurement, such that 

r -| t 

(3= b 0 0 0 

where b is an unknown constant. Thus, the measured airspeed ( v m ( t )) is expressed as 

v m = v + b (43) 

where v ( t ) is the true airspeed. The true airspeed is related to the ground speed and the wind speed as 

v = v g - v w (44) 

where v g ( t ) and v w ( t ) are the ground speed and the wind speed components along the flight path, respec- 
tively. Therefore using (44) in (43), we have 

v m = Vg - v w + b = Xp - v w + b 

where x p ( t ) is the ground position. Thus, the ground speed ( v g = x p ) is 

x p = -b + v m + v w . (45) 

Suppose an independent bias-free measurement y 0 (t) of the ground position is available (e.g., using GPS 
signals) as 

Vo = x p 

and in addition, the wind speed v w (<) measurement is available (e.g., from satellite data 15 ). Then, (45) 
can be augmented with 6 = 0 to achieve an observable system similar to (30) and the unknown bias in 
true airspeed can be estimated. Closed-loop signal boundedness and asymptotic tracking can be shown in a 
manner similar to Section V.A. 

Remark: The rate-gyro and the airspeed bias observers can be used simultaneously to estimate and 
accommodate the biases occurring simultaneously in the pitch-rate and airspeed measurements. Closed-loop 
signal boundedness and asymptotic tracking can be easily shown by extending the development in Sections 

V. A and V.B. 

VI. Application Example- Large Transport Aircraft 

Simulation studies are performed on a fourth-order longitudinal dynamics model of a large transport 
aircraft in a wings- level cruise condition with known nominal trim conditions. In the presence of unknown 
sensor bias, the control objective is to track the response of a reference model to a reference input command 
superimposed with white noise (as shown in Figure 1). As in most MRAC schemes, the adaptation/estimation 
gains (T's and rf) are selected by trial and error. 

VI. A. Dynamic Model 

The state variables for a fourth-order longitudinal dynamics model are: true airspeed v (m/s), angle-of-attack 
a (deg), pitch angle 9 (deg), and pitch rate q (deg/s). The actuators are elevator and engine throttle input, 
which produce control inputs u e (deg) and u t , respectively. The elevator input u e represents the elevator 
position (deg) and the throttle input Ut represents the thrust multiplied by a constant scale factor, and hence 
no units are used for ut . Actuator dynamics are not considered. The linear time-invariant plant is described 
by (1), where the system matrices are 

‘-0.0062 -0.0815 -0.1709 

. _ -0.0344 -0.5717 0 

0 0 0 

_ 0.0115 -1.0490 0 
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iT 


1 T 


The state is x(t) = 


v a 9 q 


, and the control input is u(t) = 


u t 


(For other flight conditions such as coordinated turns, the model would be different). 


VI. B. Reference Model 

The reference model is chosen as: 


A m — A + BK i , B m — BK 2 

where K\ is the LQR gain designed to obtain satisfactory closed-loop response. The gain I\ 2 = I 2 , so that 
B m = B. In practice, K 2 is chosen to provide proper scaling of the reference input r(t). 


VI. C. Simulation Cases 

To test the performance of various adaptive control laws presented in this paper, simulations are performed 
for the following four cases in the presence of unknown sensor bias: 


• Case 1: Standard MRAC with biased measurements (Section II) 

• Case 2: Feedback-based bias estimation with MRAC (Section III) 

• Case 3: Use of asymptotic bias estimator with MRAC (Section IV) 

• Case 4: Use of non-model-based observers with MR AC (Section V) 

In each of these cases, it is assumed that the unknown bias appears at t = 0. 

Case 1: The unknown constant bias in the state measurement is arbitrarily chosen as 


/?= 5 2 -1 10 


1 t 


(46) 


The units of the components of (3, which represent the biases in ( v,a,9,q ), are m/s, deg, deg, and deg/s 
respectively. Using the biased measurements, the standard MR AC control law (5) of Section II is implemented 
with the adaptive gains K\ and K 2 initialized to half of their true values, and the adaptation gains are chosen 
arbitrarily, so that 

K\ (0) = 0.5AA, K 2 (0) = 0.5AT 2 , F-, = 0.005/ 4 , T 2 = 0.005 1 2 . (47) 

The plant and reference model states, tracking errors, and control input are shown in Figures 2, 3, and 4, 
respectively. As stated in Section II, the closed-loop signal boundedness or bounded tracking error cannot 
be proved. However, in this example, the tracking errors appear to approach some non-zero constant values. 

Case 2: For the sensor bias as in (46), the control law in (10) is implemented where the bias estimate 
is generated using the bias estimation law in (19). The gain k 3 is initialized to zero with adaptation rate 
r 3 = 0.005/2, while other gains are chosen as in (47). The plant and reference model states, tracking errors, 
control input, and bias estimates are shown in Figures 5, 6, 7, and 8, respectively. Except for velocity, 
the controller and the bias estimator effectively compensate for the sensor bias, and achieve an acceptable 
tracking performance. The bias estimates converge to values that are not necessarily the true values. A 
better tracking performance and bias estimation may be possible if the persistent excitation (PE) condition 
is satisfied. 

Case 3: For the same sensor bias as in (46), the adaptive control law (10) of Section IV is implemented. 
The adaptive gain initialization and adaptation gains are chosen same as in Cases 1 and 2. The matrix Ap 
of the asymptotic bias estimator in (22) is arbitrarily chosen as 

Ap = 1 . 1 A , m . 

The bias estimate converges to its true value as shown in Figure 12. The plant and reference model states, 
tracking errors, and control input are shown in Figures 9, 10, and 11, respectively. The tracking errors 
converge to zero asymptotically. 

Case 4: The bias vector considered in this case is different compared to the other cases because in order 
to use the non-model-based observer developed in Section V, it is assumed that the bias exists only in the 
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airspeed and pitch-rate measurements, while reliable bias-free measurements of angle-of-attack and pitch 
angle are available. Therefore, 

r i T 


(3 


25 0 0 10 


(48) 


The adaptive gain initialization and adaptation gains are chosen same as in Case 2. The wind speed is 
assumed to be known as v w = 10 (m/s). The observer gain is chosen so as to obtain a fast convergence 
(within 1.0 sec), as indicated by the bias estimates convergence in Figure 16. The plant and reference model 
states, tracking errors, and control input are shown in Figures 13, 14, and 15, respectively. The tracking 
errors converge to zero asymptotically. 

However, if the standard MRAC control law in (5) is implemented with a bias in airspeed and pitch 
rate feedback without compensating for the sensor bias via an observer, then there is a large offset in the 
airspeed, angle-of-attack and pitch angle as indicated by the plant response shown in Figure 17. 


VI. D. Discussion 

The steady-state tracking errors improve as the design progresses from Case 1 to Case 3 (cf. Figures 3, 6, 
10). Also note the improvements in the steady-state condition of the airplane as the design progresses from 
Case 1 to Case 3 (cf. Figures 2, 5, 9). In Case 1, using just the standard MRAC control law, the airplane 
trims 5 (m/s) too slow (consistent with the 5 (m/s) too high sensor measurement) and at a flight-path 
angle (7 = 9 — a) of about —4 (deg) (descending steeper than a normal landing approach). In Case 2, with 
the addition of feedback-based bias estimation to the MRAC control law, the airplane flight-path angle is 
partially corrected to about —1 (deg) while the speed error is not affected. In Case 3, with an asymptotic 
bias estimator, both the flight-path angle and the speed are fully corrected. 

In Case 4, the bias only exists in airspeed and pitch-rate measurements. The non-model-based observer 
developed in Section V combined with MRAC is able to achieve asymptotic tracking (Figure 13). However, 
if a standard MRAC approach is used, then the steady-state response has large offsets (Figure 17). In this 
case the aircraft has a glide path angle of about —9 degrees. The aircraft is rapidly coming down while the 
reference model is not. This shows that when using standard MRAC, the sensor bias has adverse effects on 
both the steady-state tracking performance and the steady-state flight condition of the airplane. 


VII. Concluding Remarks 

This paper addressed the problem of accommodating unknown sensor bias in control of uncertain systems. 
Such a bias can have serious detrimental effects on closed-loop stability and performance. In particular, the 
standard model reference adaptive control (MRAC) law for state tracking cannot guarantee signal bound- 
edness or boundedness of the tracking error if biased sensor data is used to synthesize the control input. 
Therefore, a modified MRAC control law was developed, which incorporates sensor bias estimation along 
with gain adaptation, and was shown to guarantee signal boundedness and bounded state tracking error. 
For a case wherein a separate asymptotically convergent bias estimator is available, an MRAC control law 
was developed and was shown to provide signal boundedness and asymptotic tracking. Such asymptotically 
convergent estimators can be designed if additional independent sensor measurements are available. An 
example of a bias estimator was presented for longitudinal aircraft dynamics with constant biases in the 
pitch rate gyro and the airspeed measurements. It was also shown that the bias can be estimated and used 
in conjunction with MRAC to guarantee asymptotic state tracking and closed-loop stability. Numerical 
examples were presented to demonstrate the proposed schemes. 

The examples presented in this paper consisted of longitudinal dynamics in a wings- level cruise flight. 
Other flight conditions such as coordinated turns should be considered in future work using a six degree-of- 
freedom model. Application of the method to a full nonlinear model should also be investigated. Future work 
should also focus on extending the results to a broad class of sensor faults (including unknown scale factor 
and time- varying bias), and on methods for constructing bias estimators using data fusion for all available 
multiple redundant sensors. In addition, the methods need to be extended to address actuator faults that 
occur simultaneously with sensor faults, in order to effectively prevent loss of control in aircraft. 
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Figure 1. Reference command for all cases 
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Figure 2. Case 1: Plant and reference model states 
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Figure 3. Case 1: State tracking error 
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Figure 4. Case 1: Control input 
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Figure 5. Case 2: Plant and reference model states 
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Figure 6. Case 2: State tracking error 
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Figure 7. Case 2: Control input 
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Figure 8. Case 2: Bias estimate 
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Figure 9. Case 3: Plant and reference model states 
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Figure 10. 


Case 3: State tracking error 
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Figure 11. Case 3: Control input 
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Figure 12. Case 3: Bias estimate 
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Figure 13. Case 4: Plant and reference model states 
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Figure 14. Case 4: State tracking error 
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Figure 15. Case 4: Control input 
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Figure 16. Case 4: Bias estimate 
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Figure 17. Case 4: Plant and reference model states without using a non-model-based observer 
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